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FOREWORD 


This report presents results of an investigation of relaxation 
turbulence models. The work was done by personnel of DOW 
Industries for the Thermo- and Gas-Dynamics Division of the 
NASA Ames Research Center under Contract NAS2-8995. The 
technical monitors for the study were Mr. M. Inouye and 
Mr. R.W. MacCormackj Computational Fluid Dynamics Branch, 

NASA Ames Research Center, Moffett Field, California. 
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ABSTRACT 

Relaxation turbulence models have been intensively studied. 

The complete time-dependent mass-averaged Navler-Stokes equa- 
tions have been solved for flow into a two-dimensional compres- 
sion corner. The new numerical scheme developed by MacCormack 
has been incorporated into the developed computed code with 
an attendant order-of-magnltude reduction in computation time. 
The study consists of three interrelated parts. First, com- 
puted solutions are compared with experimental measurements 
of Law for supersonic flow for a Mach number, M^, of 2.96 and 
Reynolds number, Re^^, of 10^, with wedge angle 25°. Details 
of the relaxation process have been studied; several different 
relaxation models, including different relaxation processes 
and varying relaxation length, are tested and compared. Then 
a parametric study has been conducted in which both Reynolds 
number and wedge angle are varied. To assess effects of 
Reynolds number and wedge angle, the parametric study Includes 
the comparison of computed separation location and upstream 
extent of pressure rise ; numerical results are also compared 
with the measurements of surface pressure, skin friction and 
mean velocity field for M^=2.86, Re^^ =1.3*10® and a=24°, by 
Settles, Vas and Bogdonoff. Finally, the relaxation model is 
applied to hypersonic flow cases in which M^=8.66, Rej^=2.2»10^ 
with a highly cooled wall and wedge angles a = 27° and 36 °. 

In general, the equilibrium model predicts a sharp Increase in 
the outer-layer eddy viscosity near the separation shock; by 
contrast, the relaxation process substantially decreases the 
eddy viscosity and predicts better upstream pressure propaga- 
tion. The simple eddy viscosity model predicts, with very 
good agreement, upstream boundary conditions even in the hyper- 
sonic flow with heat transfer. However, it under-predlcts the 
eddy viscosity, and hence the skin friction, in the reattachment 
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and downstream recorapression regions. The relaxation process 
reduces the eddy viscosity further ^ and shows unfavorable 
effects on heat transfer, skin friction and reattachment 
points in these regions, thus indicating a need for further 
development and improvement of the relaxation concept. 
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NOTATION 


SYMBOL 

A 

% 

c 

V 

D 

E 


H 


h 

K 

L 


o 

M 

00 

Pr 

P 

Q 

Q ,Q 

-V 

q 

Re^ 

T 


DEFINITION 

Damping length for sublayer eddy viscosity 
Pressure coefficient, 2p^/p^u^ 

Heat transfer coefficient (Stanton number) 

Skin friction coefficient, 2x^/p^u^ 

Specific heat at constant volume 
Van Driest damping factor 
Total energy, p[e + 0 . 5 (u^ + v^ )J 
Specific internal energy, C^T 

Vector fluxes in Reynolds mass-averaged Navier- 
Stokes equation 

Vertical dimension of computational domain 
Vertical dimension of fine mesh 
Von Karman*s constant, 0.-4 

Characteristic length between leading edge and 
wedge 

Upstream boundary location 
Preestreara Mach number 
Molecular Prandtl number, 0.72 
Turbulent Prandtl number, 0.90 
Average static pressure, p = PRT 
Heat flux vector, Q^e + Q,e 

Heat flux components In Cartesian coordinates 
Velocity vector, ue^ + ve^ 

Reynolds number based on characteristic length L 
Absolute temperature, °R 
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SYMBOL 
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Ax^ 
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eq 


y 
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DEFINITION 

Time 

Vector component of conserved properties in 
Equation (1) 

Velocity at the edge of boundary layer 
Friction velocity, u^ = /x, /p„ 

Velocity components in Cartesian coordinates 
Initial relaxation location. Equation (5) 

Distance along the wall 
Cartesian coordinates 
Upstream pressure rise length 
Separation length 
Dimensionless y distance, Y"*" “ 

The y*^ of the first mesh point 
Wedge angle 

Intermittency correction factor 
Specific heat j?atio, 1.4 
Local boundary-layer thickness 

Boundary-layer thickness at the initial location x^ 
Kinematic displacement thickness 

Turbulent dynamic eddy viscosity used in calculation 
Local equilibrium eddy viscosity. Equation (4) 
Molecular viscosity 
Molecular viscosity at the wall 
Mass-averaged specific density 


vii 



H44 


N OT ATI ON (concluded) 


SYMBOL 


T 


W 


X 


DEFINITION 

Stress tensor, see Equation (1) 
Shear stress at the wall 
Relaxation length 
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1. INTRODUCTION 


Development of turbulence models is a prime factor for flow- 

field simulation. For strong nonequilibrium turbulent flows, 

superiority of current one- or two-equation models over zero- 

l-'R 

equation models remains unproven; more sophisticated models 
such as those based on the Reynolds-stress equations and sub- 
grid turbulence formulations require large (possibly imprac- 
tical) computation times. Hence, development of a modified 
simple eddy viscosity model for nonequllibriura turbulent flow 
is a very attractive approach for solving the Navier-Stokes 
equations if the concept is at all feasible. 

4 

The idea of nonequilibrium turbulence was realized by Jacobs 
in 1939 and also by Prandtl^ in 19^5- Prandtl^ studied fully 
developed turbulence and pointed out that the state of turbu- 
lence at any given moment consists partly of a residual con- 

g 

tributlon from the previous flow history. Rotta concluded 

from experimental data that when turbulent flow is perturbed 

from its equilibrium state, it takes a distance about one 

order of magnitude greater than its boundary layer thickness 

to attain a new equilibrium state. To account for this effect, 
7 

Shang and Hankey used an eddy viscosity relaxation model and 
solved the Navier-Stokes equations for a compression corner. 
Their results were sufficiently promising to encourage fur- 
ther research aimed at generalizing the simple eddy viscosity 
model to include previous turbulence history for nonequilibrium 
supersonic turbulent flows subjected to strong adverse pressure 
gradients . 

o 

In more recent applications, Baldwin and Rose, and Shang, 

9 

Hankey, and Law^ also used the relaxation model to simulate 
shock Impingement problems, and found results which were a 
significant improvement over those obtained with the simple 
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eddy viscosity model. Deiwert^^ extended the relaxation model 
to simulate a transonic flow over a thick airfoil, and Mateer, 
Brosh, and Vlegas^^ used a similar model to calculate a tran- 
sonic normal shock-wave boundary-layer Interaction. Even 
though the latter results qualitatively represent the Inter- 
active features, Mateer, Brosh, and Vlegas^^ Indicated some 
unfavorable effects on skin friction attending use of the 
relaxation model. 

The objectives of this study are: to understand how the 

relaxation model works; to study several different relaxation 
processes; and to examine the applicability of the model for 
high Mach number flows with heat transfer. The end result 
will hopefully provide some Information for further turbulence 
model development . 
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2. ANALYSIS 


2.1 G OVERN ING EQU ATI ONS 


The viscous-lnviscld Interaction Tlowrield is governed by the 
time-dependent, compressible Navier-Stokes equations. Casting 
in terms of Reynolds mass-averaged variables, and with the bulk 
viscosity and the specific turbulent energy in the normal stress 
components omitted, the resulting mean conservative relations 
are the same as their laminar flow counterpart, except for the 
addition of the Reynolds stress tensor and the Reynolds heat 
flux. Turbulent closure is accomplished by expressing the 
Reynolds stress tensor in terms of the product of eddy viscosity 
e with the mean velocity gradients. Also, a turbulent Prandtl 
number, , is used for the Reynolds heat flux. Thus, the mean 
flow equations in two dimensions can be written in integral form 
as 


where 


— f 

at 


U dVol + 


volume 


U = 



J. 


H • n dS = 0 


surface 


H = 


pq 

puq + T • e 

pvq + T • e 
Eq + T • q + Q • 
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E = p[e + 0.5(u^ + v^)] 


and e 3 e are unit vectors of the orthogonal (x,y) coordinate 

^ y ^ 

system, and n a unit normal vector of the surface enclosing the 
volume of integration. The equation of state relates the pres- 
sure, p, and density, p, to temperature, T, and specific internal 
energy, e. The perfect gas relations are 


p = pRT and e = C^T 


The molecular viscosity coefficient, y, is assumed to be a 
function of temperature only, and is evaluated by Sutherland’s 
semi-empirical formula 

y = 2.270*10“® T*^V(T+ 198.6) (slug/ft* sec ) 


In the present calculation, molecular Prandtl number, Pr, Is 
assumed to be 0.72 and turbulent Prandtl number, is taken 

to be 0.90. 


2.2 COMPUTATIONAL DOMAIN 

Figure 1 shows the computational domain and appropriate boundary 
conditions. The mesh is equally spaced in the x-dlrectlon. In 
the y-direction, a fine, geometrically-stretched mesh spacing 
is used in the region near the wall (y£h^) for resolving the 
viscous layers a coarse, equally-spaced mesh is used in the 
outer region (h£.<y<h), where viscous effects are negligible. 
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UPSTREAM 
BOUNDARY 
X CONDITIONS 



Figure 1. Computational domain and boundary conditions. 





The upstream boundary is located at a distance several boundary- 

layer thicknesses ahead of the interaction region defined from 

when the pressure starts to increase until it reaches its final 

value. The upstream conditions can be either generated from 

the present program by solving the Navier-Stokes equations over 

a flat plate that includes leading edge effects, or supplied 

from a conventional boundary- layer solution, such as Marvin and 

12 

Sheaffer’s boundary-layer code. The downstream boundary is 
positioned far enough from the corner so that all the gradients 
in the flow direction can be set to zero. Though this condi- 
tion is not exact, the boundary layer in the vicinity of this 
exit is parabolic, and the remainder of the flow is supersonic; 
hence, this condition is not expected to introduce significant 
error in the region of interest upstream. The upper boundary 
is specified by the freestream conditions. The wall surface 
is assumed impermeable, and no-slip boundary conditions are 
applied. The wall is treated either as isothermal or adia- 
batic, and the wall pressure is evaluated from an approxima- 
tion of the y-momentum equation at the wall. During the calcu- 
lation on the inner mesh, transport and stress at the Internal 
boundary, y = h^, are saved, and their net quantities are then 
used as boundary conditions for the outer mesh flowfield. 

The numerical method and special numerical procedures used in 
the calculation presented herein are described in detail in 
References 13 and l4. The computer code was initially checked 
out for laminar flow cases to ensure that the numerical scheme 
is accurate and efficient. (Details of laminar flow results 
are presented in Reference l4.) The same code is used herein 
to test the adequacy of turbulence models. During the course 
of the study the new numerical method developed by MacCormack ^ 
was incorporated in the computer code. This new method reduced 
computer times by an order of magnitude, thus facilitating 
economically feasible parametric studies (c.f.. Subsection 3*5). 
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2.3 EDDY VISCOSITY MODEL 

In the present study, a two-layer equilibrium eddy viscosity 
model ("equilibrium model") Is applied In the following 
fashion: In the Inner region, the Prandtl mixing-length 

concept Is used 

- (E)f < = ) 

where K Is the von Karman's constant (0.4) and D Is van Driest 's 
damping factor 

D = l-exp(-yZA) , A = 26 u/(/ft^Jp) 

In the law of the wall for the outer region, the modified 
Clauser's defect law is employed 

Suter = pV*^ ^3) 

where 3 is the Klebanoff intermittency correction factor 

3 = [l+ 5.5(y/6)®]"' 


and 


u 


e 


u at y = 6 (edge of boundary layer) 


6 * = 




dy 


6 = minimum of ( 6 ^ 562)5 


6^ = y|[max(pu), J < 

^2 “ u|[(Au/u) < 0.02] 


The quantity y. Is the location of zero velocity and Is a 

o max 

preassigned value, ensuring It Is located outside the boundary 
layer. Then the eddy viscosity Is determined by 
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Note that the local values of y and p are used in A, as opposed 
to the more conventional format of using the wall values ^ a 
difference which might become significant in highly cooled or 
heated wall calculations. Use of as the lower bound for 
the integration of 6* avoids picking up the unwanted negative 
(reverse) flow effects in evaluating the displacement thickness. 
The main difference of the presently employed equilibrium model 

7 

from that used by Shang and Hankey is that the velocity at the 

edge of boundary layer u^ (instead of the maximum velocity u ^ ) 

e max 

is used in the calculation of outer-layer eddy viscosity , 

^outer’ and the displacement thickness ^ 5*. Although it is dif- 
ficult in viscous-inviscid interactions to determine the edge 
of the boundary layer , care should be taken to determine u^. 
Without this effort, the maximum velocity, in general, will be 
the frees tream velocity, instead of the more typical boundary- 
layer edge velocity behind the Induced shock. Use of the free- 
stream velocity will overestimate 6*, and, in our experience, 
may cause up to a factor of two-or-more difference in the cal- 
culation of outer-layer eddy viscosity, 
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3. COMPARISON AND DISCUSSION OP RESULTS 


3.1 SUPERSONIC PLOW APPLICATIONS 

The first case computed was for M^=2.96 supersonic flow over 
a 25 ° wedge with an adiabatic wall. The flow conditions cor- 
respond identically with one of the cases studied by Shang 

rj -1 ^ 

and Hankey ' and with experimental data measured by Law. 

Flow conditions are as follows: 

= 2.96 Re^^ =10^ a = 25 "" 

= 177°R L = 1 ft adiabatic wall 

The computational domain extends from x/L= 0.795 to x/L= 1.205 
with the fine mesh boundary placed at h^= 0.025 ft and the 
outer boundary at h= 0.09 ft. In the x-dlrection 63 mesh 
points are used; in the y-direction, a total of 31 mesh points, 
20 for the fine mesh and 11 for the coarse mesh, are used. 

The mesh spacing varies from ( Ay = 0 . l48 -10" ^ft to 
Ay = 0. 446*10~^ ft for the fine mesh, and near the wall the 
mesh cells are very elongated with a ratio of Ax : (Ay)^^^ = 44 :1. 
For the case of interest, the compression is moderate (the 
pressure rise is about a factor of 5); the boundary-layer 
thickness is about the same before and after the interaction. 
Therefore, (Ay)j^j_^ is fine enough to ensure that the first 
mesh points near the wall are in the sublayer (y^ < 10) both 
upstream and downstream of the compression. 

The upstream conditions are generated with the present numerical 
code by solving the flow over a flat plate with a leading edge. 
Figure 2 shows a u-velocity profile compared with a boundary- 

1 p 7 

layer solution and Shang and Hankey 's result. The agree- 
ment Is very good. Plotted on the nondlmenslonal scale of 
(u/u^) vs y'*’, the profile Is In surprisingly good agreement 
with the law of the wall (Plgure 3). The value of y"*" of the 
first mesh point of the present study Is about 5-8 and that of 
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y X 10^ ft 



Figure 2. Plat-plate solution; M^=2.96, 

Rej^= 10^. 
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X/L = 0.795, R6l =10' 
ADIABATIC WALL 
o PRESENT RESULT 


U/U-r = 


.^^^u/UT=i7lny'*'+C 


j<— I St POINT OF SHANG a HANKEY 
^1 St POINT 




Reference 7 is Just under 10. Figures 4(a) and 4(b) show 
comparisons of pressure and skin friction calculated using 
the present equilibrium model with the experimental measure- 
ments of Law.^^ The equilibrium model tends to overpredict 
the surface pressure In the reattachment and downstream regions, 
and underpredicts the upstream pressure rise and, thus, the 
location of the separation point. Also shown In these two 
figures Is a comparison of the equilibrium results of Refer- 
ence 7 with those of the present study. Even with the good 
agreement found upstream (Figure 2), the results here are sub- 
stantially different because of the difference In formulating 
the eddy viscosity models. The use of u , Instead of u , In 
Reference 7 tends to provide a higher eddy viscosity and thus 
more mixing Inside the boundary layer, which makes It more 
resistant to separation. The modeling difference Is relatively 
Insensitive In the flat-plate solution, even though the bound- 
ary-layer profile In Reference 7 does show slightly more momen- 
tum from the Increased mixing. In the Interaction region where 
a shock wave lies above the boundary layer. It causes a 
significant difference. 

3.2 UPSTREAM RELAXATION CONCEPT 

To account for the upstream turbulence history effects, Shang 
and Hankey used a relaxation eddy viscosity model (herein 
designated "upstream relaxation" or "relaxation"), l.e., 

" = "eq - (^eq-"o)®^p(--T^) (5) 


where 

e = turbulence dynamic eddy viscosity 

e = the local equilibrium eddy viscosity evaluated 
^ from Equation (4) 

Eq = the eddy viscosity at upstream location x^ 
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(a) Surface pressure. 



L 


(b ) Skin friction. 

Figure Comparison of the results of equilibrium 

models; M =2.96, a=25°. Re, = 10’. 

oo ^ JL 
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and \ is the relaxation length. Conceptually Equation (5) 
approximates the experimental observation that. In an abrupt 
disturbance of a turbulent flow, the Reynolds stress remains 
nearly frozen at Its Initial value while It Is being convected 
along streamlines, and then exponentially approaches a new 
equilibrium state. In a numerical calculation, the Initial 
location of disturbance from which the relaxation process Is 
Initiated, x^, and a relaxation length scale which describes the 
exponential decay of the eddy viscosity distribution, X, must 
be specified. There are two limiting cases which bound the 
relaxation length. For A = 0, the turbulent eddy viscosity 
equals the local equilibrium value, and for X = », the Initial 
value Is frozen and Is used everywhere In the region x>x^. 

Using the upstream relaxation model, good Improvement In pre- 
dictions of upstream pressure rise Is obtained and a more 
pronounced pressure plateau region Is clearly demonstrated 
[Figures 5(a) and 5(b)]. For the present supersonic case, the 
relaxation length Is chosen to be X= 5<5 q, and the Initial 
location Is positioned at x^/L= 0.9387- The figures also show 
the results of equilibrium and frozen cases, which bounds the 
effect of changing the relaxation length scale from X = 0 to 
X = “. In Figure 6, the calculated density contours, obtained 
by using the relaxation model, are compared with the Interfero- 
gram and a schematic of the experimental flowfleld In Refer- 
ence 17 . Details of the Interaction pattern. Including separa- 
tion shock, separation, and reattachment shock are accurately 
simulated. 

3.3 EXPLANATION OF HOW THE RELAXATION 

MODEL WORKS 

Since the Inclusion of relaxation effects results In Improved 
comparisons of theory with experiment, It Is Important to 
understand how the relaxation model works. Figure 7 shows 
the variation of eddy viscosity at five levels from the wall. 
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(a) Surface pressure. 



Xw 

L 


(b) Skin friction. 


Figure 5. Results of different relaxation lengths; 

M^=2.96, a=25°, ReT = 10^, 6 = 0.01403 ft, 
x^L= 0.9387. ^ ° 
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SCHEMATIC DIAGRAM OF 
FLOW FIELD 



SUPERSONIC TURBULENT 
BOUNDARY LAYER OVER 
A COMPRESSION CORNER 

(Moo = 2.96, = 1Q7, a = 25°) 



INTERFEROGRAM NUMERICAL SIMULATION 

(DENSITY CONTOURS) 


Figure 6. Comparison of density contours with experimental intererogram; 

M^= 2.96, Rej^= 10^, a= 25°. 
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with the relaxation length X= 56^. The symbols represent the 

values evaluated from the local equilibrium model, 

the lines are the turbulent dynamic eddy viscosity which is 

baBed on the upstream relaxation formula. Equation (5)- The 

first three sets of data (lines and symbols) for the smaller 

values of y belong to the inner layer of the eddy viscosity 

model; the data for y/6^= 0.0928 are taken near the location 

of the transition from the inner to outer layers; and the 

outermost y-level data are calculated from the outer layer 

formula. Equation (3)- Near the wall (circle symbols. 

Figure 7) the eddy viscosity is small, e ^ < y , and is greatly 

eq w 

dominated by the wall shear stress and the local velocity grad- 
ient [see Equation (2)]. At the beginning of the interaction, 
near the wall decreases sharply because of the decrease of 
T^, reaching a minimum at separation. The minimum value of 
appears again at the point of reattachment. 

Note that = 0 at the separation and reattachment points, and 
thus, from Equation (2), D=0 and e. =0. Therefore, based 
on the local equilibrium model, [Equations (2)-(^)], there is 
no eddy viscosity, while in the real flow, turbulence fluctua- 
tions do exist at separation and reattachment. The application 
of the upstream relaxation process greatly Increases the eddy 
viscosity values in these regions and only slightly affects 
the values in the downstream region when x/L is large. 

In contrast to the effect at small y, increases sharply at 

large y because of the increase in p6* resulting from the com- 
pression of the separation shock. The relaxation process 
drastically modifies the increasing slope and provides values 
of e substantially less than that of the equilibrium model. 

The effects of the relaxation process can be examined more 
clearly in Figure 8, in which eddy viscosity profiles at three 
x-locations are plotted across the boundary layer. The reduc- 
tion of eddy viscosity from its local equilibrium value far 
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Figure 8. Eddy viscosity at three x locations. 



from the wall leads to less turbulence mixing and to less 
shear stress to balance the adverse pressure gradient, making 
the separation easier, and allowing better upstream pressure 
rise through the viscous and inviscid interaction. 

It is interesting to recall the result of the frozen case 
(X = a>) in which the initial value is used everywhere down- 
stream of x>x^. This causes the eddy viscosity in the outer 
region to be reduced by a factor of about 6 and, near the wall, 
by a factor of about 10, below the equilibrium value in the 
downstream region (see Figure 7 curves of y/6^= 0-00527 and 
y/6^ = 0. 3962 ) . Hence, while the surface pressure closely 
reaches the value of the inviscid solution, the skin friction 
is a factor of about 2 less than the equilibrium result (see 
Figure 5b ) . 

3.4 OTHER RELAXATION PROCESSES 

In addition to the upstream relaxation model, other ways of 
controlling eddy viscosity were tried. One approach is based 
on the notion that the lagging phenomena is local point-by- 
point and should be governed by the following differential 
equation 


j e - e 

d£ _ eg 

dx X 


( 6 ) 


Equation (6) is closely approximated by the following formula^^’^^ 
(herein called "point-by-polnt relaxation") 



This relaxation process avoids the question of specification 
of the upstream initial location x^, instantly accounts for 
local upstream history while the turbulence conveys downstream, 
and can handle complex cases such as multiple shock disturbances, 
etc. 
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Another observation is that a great portion of the increase 
of dynamic eddy viscosity is caused by the Increase of local 
density through shock compression. To separate the effect of 
the increase of density from the conveyance of turbulence his- 
tory, the kinematic eddy viscosity is relaxed according to the 
upstream relaxation formula 

Figures 9(a) and 9(b) compare results of these two different 
processes with the previous upstream relaxation result for 
surface pressure and skin friction. The polnt-by-point relaxa- 
tion reduces the value of equilibrium eddy viscosity more, and 
results in a larger separation region. The relaxation of kine- 
matic eddy viscosity [Equation (7)] has the opposite effect 
and is in very close agreement with the result of upstream 
relaxation in the downstream region where the relaxation 
effects are small. 

The process of relaxing only the outer layer [Equation (3)] 
has also been tested. From experimental observation, it has 
been shown that the inner layer of a boundary layer subjected 
to any disturbance adjusts quickly back to local equilibrium, 
while the outer layer may take a much longer distance to 
return to a new state of equilibrium. The process of relaxing 
only the outer layer is an extreme case of varying the relax- 
ing length across the boundary layer. It was found, however, 
that the result from relaxing only the outer layer agrees very 
closely with that from relaxing the whole boundary layer. 

3.5 PARAMETRIC STUDY AND DETAIL COMPARISON 

Instead of comparing the experimental surface pressure data 
of Law for six different wedge angles, the calculations were 
compared with newly available and more detailed data by 


21 




(a) Surface pressure. 



L 


(b) Skin friction. 


Figure 9. Results of different relaxation processes. 
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Settles, Vas and Bogdonoff. The experiments were well- 
documented, including the measurements of surface pressure, 
skin friction and detailed mean velocity field, for a wide 
range of Reynolds numbers and shock strengths. Over 50 solu- 
tions have been obtained for the entire flowfield varying 
both Reynolds number and wedge angle. The results of these 
solutions are compared with corresponding experimental data 
in Figures 10-13. Figures 10(a) and 10(b) compare the separa- 
tion locations while Figures 11(a) and 11(b) illustrate the 
upstream extent of the increased pressure level. Referring 
to the separation location comparisons. Figure 10(a) shows that 
for a wedge angle of 20° there is excellent agreement between 
the calculated values (with the equilibrium model) and experi- 
ment over the entire Reynolds number range. For the lower 
wedge angles the calculations show a larger separation distance 

than experiment. A crossplot of data at Re§ = 1.33*10^ £Fig- 

o 

ure 10(b)] shows that the calculated separation distance varies 
smoothly with wedge angle while the experimental data show an 
abrupt change near a=l6°. The calculated upstream extent of 
the increased pressure level is compared with experiment in 
Figures 11(a) and 11(b). As shown, predicted trends with 
Reynolds number agree with measurements for either model. How- 
ever, only the relaxation model predicts the correct magnitude. 

Figures 12(a) and 12(b) compare predicted and measured surface 

pressure and skin friction for the M =2.85, Re^ =1.3*10® 

^ o 

and a=2^° case. The numerical results employ the three above- 
mentioned turbulence models and a fourth model combining the 
two modifications. Details of the "pressure gradient" model 
are discussed in Reference 19. Comparing the skin-friction 
distributions [Figure 12(b)], both the equilibrium and relaxa- 
tion models substantially underpredict the skin friction down- 
stream of reattachment. The calculated reattachment locations 
are also incorrect. While the calculations employing the 
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(a) Effect of Reynolds number on Ax„ . 

o 



a, degrees 


(b) Effect of wedge angle on Ax_ . 

s 

Figure 10. Comparison of computed and measured 
separation length. Ax . 
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(a) Surface pressure. 



( ) /^ 

(b) Skin friction. w o 


Figure 12. Comparison of computed and measured surface 
pressure and skin friction for the various 
relaxation models; Mot =2.85, 5 = .83 in, 
a = 24®, = 1.3*10®. 

26 


(inch) 



>3 




(u COS a + V sin a) (ft/s ) 


(u cos a + V sin a) ( ft/s ) 


(b) Near reattachment . 


(c) Far downstream. 


Figure 13. Comparison of computed and measured mean velocity 
field. 
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pressure gradient corrections do a good job of predicting 
skin friction far downstream of reattachment, the reattachment 
location itself is still Incorrectly computed although in 
better agreement with experiment than for the equilibrium 
model. Additional comparisons with velocity profiles at three 
locations are shown in Figures 13(a)-13(c). As shown, in the 
bubble [Figure 13(a)] the experimental data are reasonably 
well predicted by the equilibrium and relaxation models , while 
the model employing only the pressure gradient correction 
underpredicts the bubble height. Far downstream at (x^“L)/6^ = 4 
[Figure 13(c)], all of the turbulence models give essentially 
the same results — all in substantial disagreement with 
experiment . 


3.6 HYPERSONIC FLOW APPLICATIONS 

At hypersonic speeds, compressibility makes drastic changes in 
the nature of turbulent flow. The turbulence disturbances are 
mostly composed of density and temperature fluctuations. At 
high hypersonic speeds, a slight local velocity change in the 
fluid would cause a large change in its temperature. Therefore, 
one might expect the test of turbulence model accuracy to be 
more severe in this case. Recall that the simple equilibrium 
eddy-viscosity model was originally derived based on the turbu- 
lent velocity fluctuations of incompressible flow. The purpose 
of this section is to investigate how the simple eddy viscosity 
model with relaxation behaves at high Mach numbers with heat 
transfer. 

The experiments selected for comparison were conducted by 
20 

Holden and measurements were made of pressure, skin friction, 
and heat transfer along the wall. The flow conditions are as 
follows : 

= 8.66 Re. = 2.2*10'^ T = 537°R 

oo J_i W 

T^ = 119.5°R L = 2.25 ft a = 27° and 36° 
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The wall is assumed to be isothermal and highly cooled; the 
ratio of wall temperature to adiabatic temperature is about 
T,/T^, =0.28. As in the previous supersonic problem, a mesh 
with 63 points in the x-direction and 31 in the y-dlrection 
is used. Now, however, the fine mesh contains 21 points for 
better resolution. The upstream boundary is located at 
x/L= 0.881 with 23 points of equal space resolution between 
the upstream boundary and the hinge line of the wedge, which 
is x/L=1.0. The mesh cells near the wall are very elongated 
with Ax/(Ay = 226 . The computation is advanced about 
95 times in the inner mesh loop for every outer mesh timestep. 

Figures l4 and 15 show the detailed comparisons of the present 
calculated results with the experimental measurements of sur- 
face pressure, heat transfer, and skin friction for ot = 27° 
and 36 °. The agreement is very good upstream of the inter- 
action region. This Implies that, at least in this study, 
the equilibrium eddy viscosity model can be used to predict 
equilibrium hypersonic turbulent flow over a flat plate. The 
fluctuations of temperature and density may be considered as 
a simple function of the fluctuation of velocity, and the 
simple eddy-viscosity relation does provide proper turbulence 
mixing in the equilibrium hypersonic, turbulent boundary layer. 


Surprisingly enough, the results also show reasonably good 
agreement in surface pressure and heat transfer in the inter- 
action region, although the peak heat transfer at the reattach- 
ment point is underpredicted. The present heat transfer 
results are significantly different from that of the laminar 
flow case studied. Instead of very low heat transfer in the 
separation region as in the laminar case, the heat transfer, 
here, drops sharply near the separation point, rises rapidly 
to a value greater than its upstream condition, gradually 
increases in the pressure plateau region, and then sharply 
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(a) Pressure coefficient. 
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(b) Heat transfer. 



(c) Skin friction coefficient. 


Figure 14. Comparisons of the results of equilibrium and 

relaxation models with experiments for M^=8.66, 
a=27°, ReL= 2.2•10^ L= 2.25 ft. 


Xq/L= 0.9762. 


6q = 0.0323 ft. 
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(c) Skin friction. 


Figure 15- 


Comparisons of the results of equilibrium and 
relaxation models with experiments; 8 . 66 ; 
a=36°, ReL= 2.2-10^, = 537°R, L=2? in. 


x^/L= 0.9392. 
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increases under strong compression. As the separation region 
becomes larger, a heat transfer plateau region forms, similar 
to that of the surface pressure. In the experimental measure- 
ments, heat transfer .does not show a sharp decrease at separa- 
tion, and, instead of the slow gradual Increase, it, indeed, 
shows a slight decrease in the pressure plateau region. It is 
still not clear whether the sharp, deep drop of heat trasnfer 
near the separation point is physically correct or is caused 
by the shortcomings of the simply eddy-viscosity model. As 
pointed out in the supersonic case (Figure 7)5 the eddy vis- 
cosity model tends to provide a very low value of eddy viscosity 
near the wall at the separation point. The underprediction of 
peak heat transfer in the reattachment region might also be 
caused by the same shortcoming of the turbulence model. 

Unfortunately, even though the points of separation and reattach- 
ment are predicted with reasonable accuracy, the skin friction 
is grossly underpredicted in the regions of reattachment and 
downstream recompression. This strongly suggests that perhaps 
the value of eddy viscosity for the present hypersonic case is 
underpredicted in these highly compressed regions. Prom the 
experience of the previous frozen case result, it may be sus- 
pected that the predicted eddy viscosity is much less than the 
real value. 

The application of the upstream relaxation process, also shown 
in Figures l4 and 15^ does show increasing upstream pressure 
propagation. However, it also shows unfavorable effects of 
decreasing skin friction and heat transfer after reattachment, 
especially in the results of the a= 27° case, in which the 
separation region is small and the initial location x^ is posi- 
tioned very close to reattachment. This may confirm our pre- 
vious argument that the simple equilibrium eddy-viscosity model 
underpredicts the eddy viscosity value in the reattachment and 
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dovmstream recompression regions, and, moreover, the relaxation 
process reduces the eddy-viscosity value still further. 

One interesting aspect of the strong compression associated with 
the hypersonic problem is the severe restriction imposed on 
(Ay)j^ln along the v/all. As shown in Figures l4(a) and l4 (b), 
the pressure rises by a factor of about 30 for the case of 
oi = 27° and about 50 for a =36°. The thickness of the boundary 
layer after compression is one order of magnitude thinner than 
that before the interaction. To ensure enough mesh resolution 
in the dovmstream highly compressed region, (Ay) ^,. ^ for the 
first mesh cell has to be very small. Figure 16 shovrs a plot 
of y+ at the first mesh point versus x. At the upstream region, 
yj is far inside the sublayer (y^^^O.SS), stays small in the 
separation region, and sharply increases right after the reattach- 
ment up to about 7. One should bear in mind that the present y 
is evaluated based on the computed skin friction. For instance, 
if these values are calculated with the experimental skin fric- 
tion, y^ vfould increase by a factor of 

|^T^^(experiment)/T^j.(calculated)J ^ » (2)'^ 

for a= 27®, vrhich may result in a yt « 10 in the downstream 
region. Hoviever, y^ is small in the separation region and the 
boundary layer in the downstream region is parabolic; hence, a 
y^ of about 10 in the dovmstream region should primarily under- 
predlct only the local skin friction and would not be expected 
to introduce significant error in the separation region. 
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Figure l6. y'*' of the first mesh points along the 

x-direction for a =27° (relaxation 
model . 
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4. SUMMARY AND CONCLUSIONS 


Relaxation turbulence models have been intensively investigated. 
The complete time-dependent mass-averaged Navier-Stokes equa- 
tions have been solved for flows over a two-dimensional com- 
pression corner. Computed supersonic solutions have been 
compared with the experimental measurements of Law and of 
Settles, Vas and Bogdonoff; also the hypersonic flows experi- 
mentally investigated by Holden have been nvimerically simulated. 
Details of the relaxation process have been studied; several 
different relaxation models, including different relaxation 
processes and varying relaxation length, have been tested and 
compared. 

In general, the equilibrium model predicts a sharp increase of 
the outer layer eddy viscosity value near the separation shock, 
vjhile the relaxation process substantially decreases the eddy 
viscosity and more accurately predicts upstream pressure propa- 
gation. The simple eddy-viscosity model predictions are in 
close agreement vrith measured upstream flow conditions in the 
hypersonic boundary layer. Hovrever, the model underpredicts 
the eddy viscosity, and hence, the skin friction in the reattach- 
ment and dovinstream recompression regions. Application of the 
relaxation process reduces the eddy-viscosity further, and shov/s 
unfavorable effects on heat transfer, skin friction, and 
reattachment point in these regions. 

Finally, the relaxation process provides a simple technique for 
including upstream turbulence history and for controlling the 
magnitude of the eddy viscosity for flov;s subjected to strong 
and abrupt disturbances such as a shock wave. Results indicate 
that the eddy viscosity should be modified in the reattachment 
recompression region for both the inner and the outer layers ; 
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the appropriate modification will yield substantially higher 
eddy viscosity values than those obtained with the equilibrium 
model. 
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